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Distorted Relation between mRNA Copy 
Number and Corresponding Major 
Histocompatibility Complex Ligand 
Density on the Cell Surface*^ 
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The major histocompatibility complex (MHC) presents 
peptides derived from degraded cellular proteins to T- 
cells and is thus crucial for triggering specific immune 
responses against viral infections or cancer. Up to now, 
there has been no evidence for a correlation between 
levels of mRNA (the "transcriptome") and the density of 
MHC-peptide complexes (the "MHC ligandome") on cells. 
Because such dependences are of intrinsic importance 
for the detailed understanding of translation efficiency 
and protein turnover and thus for systems biology in gen- 
eral and for tumor immunotherapy in practical application, 
we quantitatively analyzed the levels of mRNA and corre- 
sponding MHC ligand densities in samples of renal cell 
carcinomas and their autologous normal kidney tissues. 
Relative quantification was carried out by gene chip anal- 
ysis and by stable isotope peptide labeling, respectively. 
In comparing more than 270 pairs of gene expression and 
corresponding peptide presentation ratios, we demon- 
strate that there is no clear correlation (r = 0.32) between 
mRNA levels and corresponding MHC peptide levels in 
renal cell carcinoma. A significant number of peptides 
presented predominantly on tumor or normal tissue 
showed no or only minor changes in mRNA expression 
levels. In several cases, peptides could even be identified 
despite the virtual absence of the respective mRNA. Thus 
we conclude that a majority of epitopes from tumor-as- 
sociated antigens will not be found in approaches based 
mainly on mRNA expression studies as mRNA expression 
reflects a distorted picture of the situation on the cell 
surface as visible for T-cells. Molecular & Cellular Pro- 
teomlcs 6:102-113, 2007. 
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Major histocompatibility complex (MHC) 1 -bound peptides 
reflect a snapshot of the dynamic protein pool (1) and can 
therefore provide insights into protein turnover and the tran- 
sient cellular proteome (2-4). Class I MHC molecules present 
peptides (8-1 1 amino acids in length) (5) from degraded pro- 
teins on the cell surface to cytotoxic T-cells to allow the 
T-cells to discriminate between self and non-self (6). In addi- 
tion to its original purpose, this detection system can be used 
for the selective removal of mutated cells displaying aberrant 
MHC-peptide complexes on their cell surface (7). Therefore 
the identification of disease-related MHC-peptide complexes 
is of utmost interest for immunodiagnostics and immunother- 
apy. Huge efforts are being undertaken to identify cancer- 
specific MHC-peptide complexes applicable for tumor immu- 
notherapy (8-10). In this context, the identification of tumor- 
associated antigens (TAAs), which are predominantly 
expressed in tumor tissue, is one strategy being followed. 
Apart from differential comparisons of tumor and normal tis- 
sue proteomes (11), mRNA-based techniques are also being 
widely applied. DNA expression cloning, for example, has 
proved to be one successful method for the identification of 
TAAs- (12), but it is being increasingly complemented by large 
scale gene expression analysis using DNA microarrays (13, 14). 

The generation of MHC-peptide complexes starts in the 
cytosol where proteins are degraded by the proteasome (15). 
The source proteins for proteasomal digestion range from 
long lived proteins over short lived proteins to defective ribo- 
somal products (DRiPs) (16), which do not reach proper fold- 
ing and functional state. DRiPs are thought to be responsible 
for more than 30% of all MHC-bound peptides (1 7). Irrespec- 



1 The abbreviations used are: MHC, major histocompatibility com- 
plex; ADFP, adipophilin; dNIC, differential N-terminal isotope coding; 
DRiP, defective ribosomal product; ERAAP, endoplasmic reticulum 
aminopeptidase associated with antigen processing; HLA, human 
leukocyte antigen; NHS, nicotinoyloxysuccinimide; RCC, renal cell 
carcinoma; TAA, tumor-associated antigen; TAP, transporter associ- 
ated with antigen processing; ANXA4, annexin A4; CCNI, cyclin I; 
CD24, small cell lung carcinoma cluster 4 antigen; EHD2, EH-domain- 
containing protein 2; HMOX1, decycling heme oxygenase 1; PIGR, 
polymeric Ig receptor; PLXNB2, plexin B2; RBBP4, retinoblastoma- 
binding protein 4; TMED10, transmembrane trafficking protein 10; 
UGT1A6, UDP-glucuronosyltransferase 1 family, polypeptide A6; 
UMOD, uromodulin; ER, endoplasmic reticulum. 
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tive of their origin, a small part of peptides generated in the 
cytosol escapes further degradation and is transported via the 
transporter associated with antigen processing (TAP) into the 
endoplasmic reticulum (1 8) where they can be further trimmed 
N-terminally by proteases such as the endoplasmic reticulum 
aminopeptidase associated with antigen processing (ERAAP) 
(1 9). Finally the peptides are loaded onto MHC molecules with 
the help of several proteins, which include tapasin, calnexin, 
and calreticulin, and are exported as MHC-peptide com- 
plexes onto the cell surface. An alteration that is thought to be 
crucial for the repertoire of MHC ligands is the induction of 
immunoproteasomal subunits by interferon -y for example. 
Upon interferon y stimulation cells exchange the three pro- 
teolytic subunits of the proteasome (PSMB1 , -2, and -5) with 
subunits of different proteolytic activity (LMP2, LMP7, and 
MECL1), a hallmark for converting normal proteasomes into 
immunoproteasomes (20, 21). 

The entire MHC pathway has been investigated in depth 
often with a focus on the therapeutic potential. MHC-bound 
peptides can be used for stimulation and activation of cyto- 
toxic T-cells that, for example, have been used successfully 
for the treatment of melanoma (22, 23). The identification of 
suitable T-cell epitopes is one of the bottlenecks in the large 
scale use of this tumor immunological strategy. Utilizable 
T-cell epitopes have to be both tumor-associated and match 
the patient's MHC repertoire. Several T-cell epitopes for T- 
cell-based tumor immunotherapy have already been de- 
scribed, but the identification is laborious, and not all in vitro 
defined T-cell epitopes are potent for triggering an immune 
response also in vivo. 

Gene chip analysis is frequently used for the identification 
of TAAs (24, 25), but it has not yet been investigated whether 
an overexpression of mRNA also results in an MHC overpre- 
sentation of peptides derived from the respective protein. 
Thus, in this study MHC peptide presentation levels were 
compared with their corresponding transcription levels. We 
analyzed human clear cell renal cell carcinoma and autolo- 
gous normal tissue in a combined transcriptome and human 
leukocyte antigen (HLA) ligandome approach. HLA:peptide 
levels were quantified by mass spectrometry using stable 
isotope labels, and mRNA levels were determined by gene 
chip analysis using DNA microarrays. 

EXPERIMENTAL PROCEDURES 

Materials-The HPLC reagents, trifluoroacetic acid, acetonitrile, 
formic acid, and HPLC water, were purchased from Merck. Peptide 
modification reagent O-methylisourea hemisulfate was purchased 
from Acros Organics. l-fpH^D^ Nicotinoyioxyjsuccinimide (light or 
heavy dNIC-NHS; D represents deuterium) was synthesized as de- 
scribed elsewhere (26). 

Elution of HLA-presented Pepf/des-HLA-presented peptides were 
obtained by immune precipitation of HLA molecules from solid tissues 
using a slightly modified protocol (27) that involves the HLA-A-, -B-, 
and -C-specific antibody W6/32 coupled to protein A-Sepharose or 
CNBr-activated Sepharose (Roche Applied Science) followed by acid 
elution and subsequent ultrafiltration. 

Modification of Peptides- Modification of peptides was carried out 
as described elsewhere (26) with slight modifications. In brief, peptide 
solutions were adjusted to pH 1 1 using 1 0 m NaOH and guanidinated 
at 65 °C for 10 min using 2.5 m O-methylisourea hemisulfate solution. 



The reaction was halted by adjusting the pH to 3, and peptides were 
desalted by using Peptide Cleanup C 18 Spin Tubes (Agilent) as de- 
scribed in the manual. Peptides were loaded three times in 200-^1 
aliquots onto the Spin Tube columns. These were then washed three 
times with H 2 0, and the peptides were nicotinylated on column for 15 
min at room temperature. 500 jul of a 20 mu light or heavy dNIC-NHS 
solution in 50 mu phosphate buffer, pH = 8.5, were thus passed over 
the C 18 material with a flow rate of 33 jil/min. After washing the 
columns as described above, aminolysis of unwanted Tyr-dNIC es- 
ters was carried out by treatment with 500 ix\ of 50% hydroxylamine 
for 10 min at room temperature with a flow rate of 50 /xl/min. Subse- 
quently the columns were washed again, and the peptides were 
eluted using 4 x 50 ^l of 50% acetonitrile, 1 % formic acid. 

Mixing of Peptides-To mix tumor and normal MHC peptides in a 
total peptide ratio of 1 :1 , absorption of isotope-labeled peptide pools 
was determined at 260 nm (28). Peptide content of the tumor and 
normal sample was calculated using the molar extinction coefficient 
for nicotinic acid in 50% acetonitrile ([cepsilon] = 1430 m _1 ), and 
equal amounts of peptides were mixed. 

Microcapillary LCMS— Peptide analysis was carried out as de- 
scribed elsewhere (26) using an Ultimate HPLC system (Dionex) with 
a gradient ranging from 15 to 55% solvent B within 170 min. Mixed 
tumor and normal samples were recorded in an LCMS experiment 
without fragmentation using a hybrid quadrupole orthogonal acceler- 
ation time of flight MS/MS system (Q-TOF, Micromass) equipped with 
a micro-ESI source. For sequence analysis tumor and normal sam- 
ples were analyzed separately in individual LCMS/MS experiments. 

Peptide Sequence Analysis and Peptide Quantification -Peptide 
sequence analysis was carried out manually. Therefore MS/MS spec- 
tra were smoothed using MassLynx 4.0 software (Savitzky-Golay, 
three smooth windows, two smooths). For manual peptide identifica- 
tion, sequence tag searches were done using Mascot 2.0 software 
(peptide and MS/MS tolerance, 0.2 Da; National Center for Biotech- 
nology Information non-redundant (NCBInr) database updated 
monthly and restricted for search to human entries; 138,263 entries 
for human proteins at the time when the searches were done), and 
relevant hits were assessed manually and not by evaluation of the 
Mascot score. Criteria for manual identification were: a reasonable 
interpretation of at least 95% of all fragment peaks, complete se- 
quence coverage with MS/MS fragments, and signal intensities of 
fragment ion peaks that match breakage probabilities of the respec- 
tive sequence. For more than 50% of all sequenced peptides MS/MS 
spectra with both light and heavy dNIC isotope label were available, 
facilitating manual sequencing of the peptides. Ions containing the N 
terminus showed a mass shift of 4 Da due to the dNIC isotope label 
(see Fig. 2, C and D). In total 363 peptides were sequenced in the 
three analyzed RCC samples. All identified peptides were blasted 
against the NCBInr database (see above), and all peptides that could 
not be unambiguously allocated to only one protein sequence were 
excluded from further analysis unless the mRNA expression ratios for 
these potential source proteins were equal. For all peptides that could 
be linked to an mRNA expression ratio Supplemental Table 1 lists 
protein information, HLA presentation ratio, mRNA expression ratio, 
HLA restriction, source, precursor mass, and charge state. All rele- 
vant MS/MS spectra used for peptide identification are available upon 
request. For peptide quantification mixed tumor and normal tissue 
samples were analyzed in an MS experiment. Peptide pairs were 
identified manually. Afterward mass chromatograms for the light and 
heavy versions of the peptide were calculated, and scans contain- 
ing the peptide of interest were summed up. For background sub- 
traction an equal number of scans with a similar retention time was 
subtracted. The relative amount of light and heavy dNIC peptides 
was calculated using the peak heights of the first isotope peaks 
(see Fig. 2B). 
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Fig. 1 . Strategy for differential 
mRNA and MHC ligand analysis. A, 

mRNA and MHC ligands were isolated 
from tumor and autologous normal tis- 
sue. Relative quantification was carried 
out for both specimens, and the correla- 
tion between the quantitative data was 
determined. B, modification strategy for 
MHC ligands. Tumor-derived ligands 
were modified with heavy dNIC reagent; 
normal tissue was modified with the light 
form. Sequence analysis was carried out 
individually for tumor and normal sam- 
ples (K* indicates a guanidinated lysine 
side chain); for ratio analysis both sam- 
ples were mixed in a total peptide ratio 
of 1 :1 . C, dNIC modification reagent with 
X representing the site of stable isotope 
labeling. D, deuterium. Depiction of kid- 
ney: inspired by Prof. Dr. K. U. Benner, 
University of Munich. 
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Assessment of False Positive Rate of Peptide Identification-To 
estimate the false positive rate of the identified peptides, a database 
was designed that contained both the European Bioinformatics Insti- 
tute International Protein Index human database (IPI, version 3.21, 
containing 60,822 entries) and the reversed European Bioinformatics 
Institute International Protein Index human database (rlPI). 50 ran- 
domly chosen peptides from RCC099 were searched using the man- 
ually identified sequence tags in this database with the Mascot soft- 
ware. In summary 20% of all searches returned only one peptide with 
a significant Mascot score, and 66% of all searches returned only one 
hit (all in IPI). Only 6% of all searches resulted in more than three 
peptide hits (both in IPI and rlPI, Supplemental Table 2). Eight of 
theses peptides and eight additional peptides, also identified in 
RCC099, were chemically synthesized and modified (Supplemental 
Table 3). Fragmentation spectra of these synthetic peptides were 
compared with the fragmentation spectra recorded from RCC099. All 
synthetic peptides showed exactly the same fragmentation pattern as 
the peptides identified in RCC099. Thus we conclude that the false 



positive rate in our sequence analysis is below 1 %. Comparisons of 
MS/MS spectra obtained from synthetic peptides and from RCC099 
are available upon request. 

Gene Expression Analysis by High Density Oligonucleotide Mi- 
croarrays—RNA isolation from tumor and autologous normal kidney 
specimens and gene expression analysis by Affymetrix Human Ge- 
nome U133 Plus 2.0 oligonucleotide microarrays (Affymetrix, Santa 
Clara, CA) were performed as described previously (29). Data were 
analyzed with the GeneChip Operating Software (Affymetrix). Pairwise 
comparisons between tumor and autologous normal kidney were 
calculated using the respective normal array as base line. 

RESULTS 

Principles of Quantitative Peptide and mRNA Analysis— The 
general strategy used to assess the correlation between 
mRNA ratios and their corresponding ratios of HLA-presented 
peptides is depicted in Fig. 1 A mRNA was isolated from three 
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Fig. 2. Peptide ratio analysis and sequence identification. A, for estimation of methodical deviations a mixture of six synthetic peptides 
was modified with the dNIC strategy using either heavy or light dNIC in two independent experiments. Heavy and light samples were mixed 
in a total peptide ratio of 0.8:1 and 1 :1 . Actual peptide ratios were determined by mass spectrometry; mean and S.D. are stated. B-D, MHC 
ligands were isolated from RCC1 00 tumor and normal tissue and subsequently labeled differentially with heavy or light dNIC, respectively. Ratio 
determination of differentially labeled ALASHLIEA was done in MS mode (B). MS/MS fragment spectra (C and D) were used for sequence 
identification. 



different RCC tumor tissues and respective autologous nor- 
mal tissues for subsequent gene chip analysis. HLA-pre- 
sented peptides were isolated from the same tumor and nor- 
mal tissue samples (30) and quantified relatively as shown in 
Fig. 1B using the differential N-terminal isotope coding (dNIC) 
strategy (26, 28). The HLA ligandome isolated from tumor 
tissue was derivatized with "heavy" nicotinic acid (deuterium- 
bearing, Fig. 1C), whereas for the HLA ligandome of normal 
tissue the "light" nicotinic acid (hydrogen-bearing, Fig. 1C) 
was used. As a consequence, all peptides from tumor tissue 
possessed the very same physicochemical features as their 
counterparts with identical sequence from normal tissue but 
differed in mass by 4 Da. To achieve the highest possible 
accuracy and sensitivity both for quantification and sequence 
identification, quantification and sequence identification ex- 
periments were carried out separately. Therefore, after deri- 
vatization one aliquot of each peptide pool was used for 
relative quantification. Heavy dNIC-peptides from tumor sam- 
ples and light dNIC-peptides from normal tissue were mixed 
in a total peptide ratio of 1 :1 and subjected to on-line coupled 
ESI-LCMS. Pairs of identical peptides derived from tumor and 
normal tissue were identified on the basis of a 4-Da mass 
difference introduced by the isotope label and due to their 
identical retention times in chromatography. For differential 
quantification, the ratio of signal intensities between monoiso- 



topic peaks of the tumor- and normal tissue-derived peptide 
was calculated (Fig. IB, right panel). For peptide sequence 
analysis, the remaining aliquot of derivatized HLA ligand pools 
was used for peptide sequence analysis by on-line coupled 
ESI-LCMS/MS. Fragmentation experiments were carried out 
separately for each tumor and normal tissue (Fig. 16, bottom 
panel). Thus an unambiguous allocation to tumor or normal 
tissue was possible for each peptide signal even in the case of 
peptides for which the corresponding sequence could not be 
determined. The separate analysis of tumor and normal tissue 
also allowed the exclusion of singlet signals, which occurred 
due to incomplete modification with dNIC, from further anal- 
ysis. Quantitative information was linked to sequence infor- 
mation by peptide individual features combining mass, charge 
state, and retention time. To estimate the intrinsic method- 
ological error, this quantification strategy was also applied to 
synthetic peptide pools, resulting in a calculated error for 
peptide pairs with high MS signal intensity of less than 10% 
(Fig. 2A). 

mRNA Versus HLA Peptide: Quantitative Comparisons-\n 
our study, this relative mRNA and HLA peptide quantification 
strategy was extended to three RCC samples sharing two 
HLA alleles: RCC099 (HLA-A*02, -A*03, -B*27, and -B*57), 
RCC100 (HLA-A*02, -A*03, -B*07, and -B*18), and RCC110 
(HLA-A*02, -A*68, -B*18, and -B*27). Table I gives an over- 
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Table I 

Absolute numbers of evaluated mRNA ratios, HLA ligand ratios, and 
peptide sequences identified in this study 



In the last column for all RCCs the number of successful allocations 
of peptide sequences to the corresponding HLA ligand ratios and 
mRNA ratios is given. 



Source 


mRNA 
ratios 


HLA 
ligand 
ratios 


Identified 
peptide 
sequences 


Allocated 
sequence, 
peptide, and 
mRNA ratio 


RCC099 


18,976 


452 


213 


166 


RCC100 


24,485 


142 


59 


43 


RCC110 


26,571 


134 


91 


64 


Total 


70,032 


728 


363 


273 



view of the identified mRNA and HLA peptide ratios, the 
number of identified sequences, and the peptides for which 
sequence information could be linked to its corresponding 
presentation ratios. In total, intensity ratios between tumor 
and normal tissue were determined for 728 peptide pairs (Fig. 
2B), 363 peptides of which could be sequenced (Fig. 2, C and 
D). 273 peptide sequences could be linked to a specific 
mRNA allowing the evaluation of a possible correlation be- 
tween differences in gene expression and HLA peptide pres- 
entation (Supplemental Table 1). 

To be able to combine the HLA-mRNA comparisons of all 
three RCCs datasets, relative mRNA expression ratios and 
HLA peptide presentation levels had to be normalized individ- 
ually for each analyzed tissue pair (Fig. 3A). To do so, histo- 
grams of iogarithmized HLA peptide presentation levels were 
plotted and subsequently fitted assuming a Gaussian distri- 
bution. Each fitted HLA histogram plot had its maximum near 
0, representing a 1 :1 ratio of presentation on tumor and nor- 
mal tissue. To take experimental inaccuracies into account 
regarding the mixing of HLA peptides from each tumor and 
normal tissue pair, the HLA presentation ratios were normal- 
ized such that the maxima of their fitted histogram plots were 
located at 0. Normalized HLA peptide presentation data of all 
three RCCs were cumulated and fitted assuming a Gaussian 
distribution to determine the presentation level representing 
the 5% highest over- or underpresentation (x T op5%> Fig. 36). 
An analogue normalization and assessment of the x TOP 5 o /o 
was carried out for the mRNA expression data. 

Finally normalized mRNA ratios were plotted against their 
corresponding normalized HLA presentation ratio (Fig. 4), and 
the Spearman rank correlation coefficient was calculated at 
r = 0.32, showing an even weaker correlation of transcrip- 
tome to HLA ligandome than the correlation of transcriptome 
to proteome (r = 0.45) (4). For qualitative data analysis and 
simplified data discussion Fig. 4 was divided into nine areas 
by the x TOP5% margins of mRNA and peptide ratios as indi- 
cated by the clashed lines (Fig. 3S). 

HLA Ligand Pools Derived from Tumor and Normal Tissue- 
For -75% of all identified HLA ligands we could not detect 
significant changes either in their HLA presentation ratio or in 




2log(-jL)mRNA ». 2ksg^j ) MHC ligand — » 

B RCC099 + RCC100 + RCC1 1 0 




2log<^)mRNA ► 2 log<£) MHC ligand— ► 



Fig. 3. Normalization of quantitative mRNA and MHC ligand 
data. mRNA ratios {left panels) and MHC ligand ratios (right panels) 
between tumor and normal tissue were logarithmized, and frequency 
counts were made with bin sizes from 0.3 to 0.7. The bin sizes were 
chosen so that the Gaussian fit resulted in the best Ft 2 values. Using 
these Gaussian fit functions, the mean value (x c ) of each data set was 
calculated. A shows data sorted for RCC099, RCC1 00, and RCC1 1 0. 
All data were normalized to x c = 0. B, the margins of relative gene 
transcript expression representing the 5% strongest over/underex- 
pression (x TOP s% ) were calculated as follows. Normalized mRNA data 
for all RCCs were cumulated and fitted assuming a Gaussian distri- 
bution, the area under this curve was calculated, and the x values 
corresponding to 5% (for underpresentation) and 95% (for overpre- 
sentation) of the area were determined. x TOP s% of the HLA ligandome 
ratios was calculated similarly. T, tumor; N, normal. 

their corresponding mRNA expression levels (Fig. 4, area V). 
Regarding mRNA expression this might have been expected 
as it indicates on a transcriptional basis that the tumor tissue 
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Fig. 4. Scatter plot representation of 
the correlation between relative 
mRNA expression and relative MHC 
ligand presentation. The log 2 of the 
mRNA and MHC ligand ratios of tumor 
versus normal tissue were calculated 
and plotted for 273 pairs of mRNA and 
MHC ligand ratios. Linear regression is 
shown as a solid line, and the Spearman 
rank correlation coefficient was calcu- 
lated. Dashed lines designate the mar- 
gins for the 5% of mRNA (x = ±1 .3) and 
peptide (y = ±1.5) specimens with the 
highest under- and overexpression, re- 
spectively, presentation. Thus, the plot is 
divided into nine areas (see Inset). Se- 
lected loci are highlighted. T, tumor; N, 
normal. 




still resembles its corresponding normal tissue. As to the HLA 
level this is more surprising as gene chip analysis showed 
up-regulation of interferon y-inducible proteins in the antigen- 
processing machinery (Table II). In all tumor tissues analyzed, 
immunoproteasomal subunits were clearly overexpressed (up 
to 11 -fold for LMP2 in RCC110). Although immunoproteaso- 
mal activity is thought to change the HLA-bound peptide rep- 
ertoire (31) and although it has been shown that both normal 
and immunoproteasomes are capable of destroying each oth- 
er's epitopes (32, 33), we could not detect major differences in 
the HLA ligandome of tumor versus normal tissue. 

Multiple HLA Ligands from One Source Profe/n-Table III 
lists mRNA and HLA presentation ratios of proteins that are 
highlighted in Fig. 4 and discussed below. Analysis of individ- 
ual peptide pairs indicated tumor-associated changes in HLA 
presentation ratios (Table IV). For 20% of all proteins with 
several identified HLA ligands, HLA presentation ratios dif- 
fered between the individual peptides more than 4-fold. For 
four peptides, N-terminally trimmed forms were identified 
(CCNI, HLA-B, RBBP4, and TMED1 0). The presentation ratios 
of these peptide variants changed up to 5-fold. This might 
indicate tumor-associated changes in antigen processing (34) 
and especially alterations in trimming by the ERAAP (35). 

HLA Ligands Identified Only in One Tissue Specimen- 
Several HLA ligands were found to be presented exclusively 
on one tissue specimen (Table V), although mRNA ratios 
changed only marginally in most cases. Peptides derived from 
UDP-glucuronosyltransferase 1 family, polypeptide A6 
(UGT1A6), uromodulin (UMOD), and the polymeric Ig receptor 
(PIGR) were presented uniquely on normal tissue (Fig. 4, areas 
VIHX); peptides derived from annexin A4 (ANXA4), plexin B2 



(PLXNB2), decyciing heme oxygenase 1 (HMOX1), and adi- 
pophilin (ADFP) as well as from the small cell lung carcinoma 
cluster 4 antigen CD24 were presented solely on tumor tissue 
(Fig. 4, areas II and III). Peptides presented specifically on one 
tissue specimen might reflect both differences in antigen 
processing between one tumor/normal tissue pair and pa- 
tient-specific variations: identical peptides from ANXA4 and 
ADFP were identified in different tumors but without being 
presented uniquely. For 60% of all source proteins in the 
group of the top 5% overpresented peptides (Fig. 4, areas 
l-llf), tumor association has been reported, and this number 
actually increases to nearly 80% for peptides presented ex- 
clusively on tumor tissue (Supplemental Table 4). As far as 
tumor immunotherapy is concerned, these uniquely pre- 
sented peptides are excellent targets, but they also underline 
the importance of patient-specific adaptations in such an 
approach (9). 

HLA Ligands without Corresponding mRNA— The fact that 
there are differences between the transcriptome and HLA 
ligandome becomes even more obvious when one considers 
that several peptides could be identified for which no mRNA 
was detectable in the respective tissues (Table VI). In 
RCC1 00, for example, the mRNA for the EH-domain-contain- 
ing protein 2 (EHD2) was detectable neither in tumor nor in 
normal tissue. Nevertheless the EHD2-derived peptide 
ALASHLIEA was identified on tumor as well as on normal 
tissue (Fig. 2, B-D). These HLA ligands are presumably de- 
rived from long lived proteins from which the mRNA has 
already been degraded (36) or represent mutated peptides for 
which the correct allocation to the correct mRNA was not 
possible any more. 
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Gene chip analysis of gene transcripts involved in antigen processing 
Analysis comprised constitutive proteasomal subunits (PSMB1, -2, and -5), immunoproteasomal subunits (LMP7, -2, and MECL1), subunits 
of the immunoproteasomal regulator PA28 (PSME1-3), proteins of the HLA peptide-loading complex (TAP1 and -2, tapasin, calreticulin, 
calnexin, and PDIA3), HLA classes I and II (HLA-A, -B, -C, -DR, -DP, and -DQ), the ERAAP, and proteins of the SEC61 translocon that can be 
involved in peptide export from the ER. Proteins inducible by interferon y are highlighted in bold, histocomp., histocompatibility; cass., 
cassette; aminopep. assoc. w. antigen process., aminopeptidase associated with antigen processing. 



Gene symbol 


Gene title 










RCC099 


RCC100 


RCC110 


PSMB1 


Proteasome subunit 01 


0.4 


0.4 


0.6 


PSMB2 


Proteasome subunit /32 


0.3 


0.1 


0.5 


PSMB5 


Proteasome subunit 03 


0.0 


-0.2 


0^2 


LMP7 


Proteasome subunit 08 


2.5 


2.1 




LMP2 


Proteasome subunit 09 


1.8 


1.7 




MECL1 


Proteasome subunit 010 
Proteasome activator subunit 1 


1.3 


0.8 




PSME1 


0.2 


0.2 




PSME2 


Proteasome activator subunit 2 


0.1 


-0.3 


1.4 


PSME3 


Proteasome activator subunit 3 


-0.1 


-0.9 


-0.6 


CALR 


Calreticulin 


0.0 


0.0 


0.1 


CANX 


Calnexin 


1.5 


0.8 


2.0 


TAP1 


Transporter 1, ATP-binding cass. 


1.5 


1.8 


2.6 


TAP2 


Transporter 2, ATP-binding cass. 


0.3 


1.1 


2.6 


TAPBP 


Tapasin 


1.4 


1.1 


1.6 


HLA-A 


Major histocomp. complex I A 


0.8 


0.3 


1.1 


HLA-B 


Major histocomp. complex I B 


1.3 


0.8 


2.1 


HLA-C 


Major histocomp. complex I C 


0.9 


0.5 


1.6 


HLA-DP 


Major histocomp. complex II DP 


1.4 


0.6 


2.5 


HLA-DQ 


Major histocomp. complex II DQ 


2.0 


1.3 


3.5 


HLA-DP. 


Major histocomp. complex II DR 


1.1 


0.6 


2.4 


PDIA3 


Protein-disulfide isomerase A3 


0.2 


0.1 


0.7 


SEC61A1 


Sec61 a1 subunit 


0.6 


-0.5 


0.5 


SEC61A2 


Sec61 <*1 subunit 


0.0 


-0.9 


-1.2 


SEC61B 


Sec61 a1 subunit 


-0.1 


-0.2 


0.6 


SEC61G 


Sec61 a1 subunit 


1.9 


0.9 


1.2 


ERAP1 


ER aminopep. assoc. w. antigen process. 


0.7 


0.7 


0.8 



Table III 

Diversity of peptide and RNA ratios highlighted by 1 1 examples 
Genes were sorted according to the area classification introduced in Fig. 4. HLA ligand ratios and mRNA ratios were calculated between 
tumor and autologous normal tissue. Gene symbol and title can be found at www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=gene. 





PLXNB2 


Plexin B2 


TYTDRVFFL 


9.58 


-0.50 


RCC110 




ADFP 


Adipose differentiation-related protein 


VRLGSLSTK 


4.22 


-0.08 


RCC099 




CD24 


CD24 antigen 


RAMVARLGL 


5.32 


0.71 


RCC100 




ANXA4 


Annexin A4 


DEVKFLTV 


3.52 


1.41 


RCC100 




ADFP 


Adipose differentiation-related protein 


SLLTSSKGQLQK 


2.32 


1.71 


RCC100 




ANXA4 


Annexin A4 


DEVKFLTV 


1.95 


2.20 


RCC110 




HMOX1 


Heme oxygenase (decycling) 1 


KIAQKALDL 


2.30 


3.70 


RCC110 


VI 


ADFP 


Adipose differentiation-related protein 


VRLGSLSTK 


-0.81 


3.60 


RCC110 


VII 


UMOD 


Uromodulin 


RAFSSLGLLK 


-5.49 


-7.19 


RCC100 


VIII 


UGT1A6 


UDP-glucuronosyltransferase 1 A6 


ALGKIPQTV 


-4.55 


1.12 


RCC099 


IX 


PIGFt 


Polymeric immunoglobulin receptor 


FSWINQLR 


-4.73 


1.72 


RCC099 



HLA Ligands Identified in Several Patient Samples- Pep- 
tides identified in more than one patient allowed the evalua- 
tion of patient-specific variations in mRNA expression ratios 
and HLA presentation ratios (Supplemental Table 5). More 
than 12% of the HLA-presented peptides exhibited a change 
in the presentation level in different patients that was greater 
than 4-fold and included many of the significantly over- or 
underpresented peptides. The peptide VRLGSLSTK derived 



from the tumor-associated antigen ADFP (37) was underpre- 
sented on RCC110, but mRNA levels were clearly up-regu- 
lated in the tumor tissue (Fig. 4, area VI). Yet the same peptide 
was presented exclusively on RCC099, although the mRNA 
levels of ADFP remained unchanged between tumor and nor- 
mal tissue (Fig. 4, area II). Another different ADFP-derived 
peptide (SLLTSSKGQLQK) was presented on RCC100 by 
HLA-A*03 with an mRNA ratio that matched its HLA presen- 
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Table IV 

Multiple peptides from one source protein identified in one tumor and normal tissue pair 
HLA presentation ratios were calculated between a corresponding tumor and normal tissue pair. Peptides that were identified also in 
N-terminally trimmed forms are marked in bold. Gene symbol and title can be found at www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=gene. 
elongat., elongation; histocomp., histocompatibility; Heterogen., heterogeneous; ribonucleoprot., ribonucleoprotein; bind, prot, binding 
protein; activ., activator. 



Gene symbol 


Gene title 


Sequence 


log 2 HLA 
ligand ratio 


log 2 mRNA 
ratio 


restriction Sourca 


ALDOA 


Aldolase A 


ALSDHHIYL 


3 38 




A02 


RCC099 


ALDOA 


Aldolase A 


RTVPPAVTGITF 




9 79 


B57 


RCC099 


CCNI 


Cyclin I 


LLDRFLATV 


0 93 


n -is 


A02 


RCC099 


CCNI 




SLLDRFLATV 




n i* 


A02 


RCC099 


EEF2 


Eukaryotic translation elongat. factor 2 


ILTDITKGV 




o 49 


A02 


RCC099 


EEF2 


Eukaryotic translation elongat. factor 2 


RRWLPAGDAL 


037 


049 


B27 


RCC099 


FLNA 


Filamin A 


GTHKVTVLF 


006 


1 09 


B57 


RCC099 


FLNA 


Filamin A 


GTHTVSVKY 


1 01 




B57 


RCC099 


FLNA 


Filamin A 


GVHTVHVTF 


-0 24 




B57 


RCC099 
RCC099 


GNB2L1 


Guanine nucleotide-bind. prot. b2-like 1 


KTIKLWNTL 


1 06 


1 02 


A02 


GNB2L1 


Guanine nucleotide-bind. prot. b2-like 1 


YTDNLVRVW 




1 09 


B57 
A02 


RCC099 


HLA-A or -G 


Major histocomp. complex I A or G 


VMAPRTLLL 


-0 23 


1 60 


RCC110 


HLA-A 


Major histocomp. complex I A 


VMAPRTLVL 


040 


1 60 


A02 


RCC110 


HLA-B 


Major histocomp. complex I B 


AAQITQRKW 


1 84 


1 29 


B57 


RCC099 


HLA-B 


Major histocomp. complex I B 


TAAQITQRKW 


161 


1 29 


B57 


RCC099 


HNRPM 


Heterogen. nuclear ribonucleoprot. M 


KSRGIGTVTF 


-0 74 


047 


B57 


RCC099 
RCC099 


HNRPM 


Heterogen. nuclear ribonucleoprot. M 


LLFDRPMHV 


043 


047 


A02 


LMNA 


Lamin A/C 


KAGQWTIW 




0 07 


B57 


RCC099 
RCC099 


LMNA 


Lamin A/C 


MRARMQQQL 




0 07 


B27 


PPP1CA or -B or -C Protein phosphatase 1a or b or c 


KYPENFFLL 


620 


0 <?5 


C 


RCC110 


PPP1CA 


Protein phosphatase 1 a 


SIIGRLLEV 


047 




A02 


RCC110 
RCC099 


RBBP4 or -7 


Retinoblastoma-binding protein 4 or 7 


HTAKISDFSW 


-0 13 


* 


B57 


RBBP4 or -7 


Ftetinoblastoma-binding protein 4 or 7 


TAKISDFSW 


043 




B57 


RCC099 


RPS16 


Ribosomal protein S16 


ISKALVAYY 


19 


0 89 


B57 


RCC099 
RCC099 


RPS16 


Ribosomal protein S16 


KLLEPVLLL 


1 rw 


nfl9 


A02 


SCO 


Stearoyl-CoA desaturase 


ARLPLRLFL 


3 96 


1 52 


B27 


RCC099 


SCD 


Stearoyl-CoA desaturase 


ITAGAHRLW 


1.50 


1.52 


B57 


RCC099 


SSR1 


Signal sequence receptor 


VLFRGGPRGLLAVA 


-2.00 


0.70 


A02 


RCC110 


SSR1 


Signal sequence receptor 


VLFRGGPRGSLAVA 


-2.42 


0.70 


A02 


RCC110 


SSR1 


Signal sequence receptor 


VLFRGGPRGLLAVA 


0.89 


0.37 


A02 


RCC099 


SSR1 
STAT3 


Signal sequence receptor 


VLFRGGPRGSLAVA 


-1.35 


0.37 


A02 


RCC099 


Signal transducer and activ. of transcript. 3 


EERIVELF 


0.10 


-0.80 


B18 


RCC110 


STAT3 


Signal transducer and activ. of transcript. 3 


EELQQKVSY 


0.38 


-0.80 


B18 


RCC110 
RCC099 


TMED10 


Transmembrane trafficking protein 10 


FLLGPRLVLA 


-0.91 


-0.08 


A02 


TMED10 


Transmembrane trafficking protein 10 


LLGPRLVLA 


1.49 


-0.08 


A02 


RCC099 
RCC099 


TMEM66 


Transmembrane protein 66 


KGWDGYDVQW 


0.39 


0.72 


B57 


TMEM66 


Transmembrane protein 66 


RRLDPIPQL 


-3.55 


0.72 


B27 


RCC099 


TNS1 


Tensin 1 


HAKVLEFGW 


0.81 


0.32 


B57 


RCC099 


TNS1 


Tensin 1 


FLIETGPRGV 


0.45 


0.32 


A02 


RCC099 


VIM 


Vimentin 


NLAEDIMRL 


0.27 


1.82 


A02 


RCC099 


VIM 


Vimentin 


NYIDKVRFL 


1.02 


1.82 


C 


RCC099 



tation ratio (Fig. 4, area III). HLA peptide ratios exhibited 
stronger patient individual changes than mRNA ratios might 
indicate. This is due to the fact that changes in the HLA 
ligandome reflect not only mRNA ratio alterations but also 
changes in the degradome, which to some extent is con- 
served in the HLA ligandome (38). 

DISCUSSION 

Stable Isotope Labeling of ex Vivo Isolated MHC Molecules- 
Stable isotope labeling is the method of choice for protein 



quantification and is widely used in proteomics. For ex vivo 
prepared MHC peptide samples the standard isotope labeling 
methods are either not feasible or are not applicable to the 
majority of MHC-bound peptides due to their amino acid com- 
position. To overcome these impediments we have established 
the dNIC strategy; to our knowledge this is the first study using 
an isotope-based quantitative MS method for a large scale ex 
vivo comparison of the HLA ligandome of two tissue specimens. 
Patient samples were chosen for analysis according to two 
criteria: categorization of the RCCs into the clear cell type by 
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HLA ligands exclusively identified in either tumor or normal tissue 
Values were calculated against background and sorted according to their HLA presentation ratio between tumor and autologous normal 
tissue. Gene symbol and title can be found at www.ncbi.nlm.nih.gov/entrez/query .fcgi?db=gene. histocomp., histocompatibility; bind. prot. 
binding protein; N/A, not assigned. 



Gene symbol 


Gene title 


Sequence 


log 2 HLA 
ligand ratio 


log 2 mRNA 
ratio 


HLA 

restriction 


ource 


PLXNB2 


Plexin B2 


TYTDRVFFL 


9 58 


0 00 




RCG1 1 0 


PPP1CA or -B or -C 


Protein phosphatase 1 a or b or c 


KYPENFFLL 


620 






RCC110 


CD24 
SLC17A3 


CD24 antigen 


RAMVARLGL 


532 


(171 


N/A 


RCC1 00 


Solute carrier family 17, member 3 


ARYGIALVL 


486 


0 62 


R97 


RCC099 


ADFP 


Adipose differentiation-related protein 


VRLGSLSTK 


4 22 


-008 




RCC099 


MAT1A or -2A 


Methionine adenosyltransferase la or lia 


RRVLVQVSY 


4 08 


~050 


B27 


RCC1 10 


SCO 


Stearoyl-CoA desaturase 


ARLPLRLFL 


3.96 


1.52 


B27 


RCC099 


ANXA4 


Annexin A4 


DEVKFLTV 


3.52 


1.41 


B18 


RCC100 


P0LR2C 


Polymerase II polypeptide C, 33kDa 


KLSDLQTQL 


3.52 


0.10 


A02 


RCC110 


HM0X1 


Heme oxygenase 1 


KIAQKALDL 


2.30 


3.70 


A02 


RCC110 


CXCL14 


Chemokine ligand 14 


RLLAAALLL 


2.09 


0.47 


A02 


RCC099 


0GG1 


8-Oxoguanine DNA glycosylase 


VLADQVWTL 


-3.07 


0.17 


A02 


RCC099 


MYL6 


Myosin light polypeptide 6 


FVRHILSG 


-3.45 


0.12 


N/A 


RCC099 


GNBS 


Guanine nucleotide-bind. prot. b5 


ILFGHENRV 


-3.72 


-0.78 


A02 


RCC099 


HLA-B or -C 


Major histocomp. complex I B or C 


DTAAQITQR 


-4.41 


1.58 


A68 


RCC110 


RPL15 


Ribosomai protein L15 


EVILIDPFHK 


-4.51 


-0.75 


A68 


RCC110 


UGT1A6 


UDP-glucuronosyltransferase 1 A6 


ALGKIPQTV 


-4.55 


1.12 


A02 


RCC099 


PIGR 


Polymeric immunoglobulin receptor 


FSWINQLR 


-4.73 


1.72 


A03 


RCC099 


UMOD 


Uromodulin 


RAFSSLGLLK 


-5.49 


-7.19 


A03 


RCC100 


MY01C 


Myosin IC 


FLDHVRTSF 


-5.69 


0.36 


A03 


RCC100 



Table VI 

HLA ligands identified from RCC099, RCC100, and RCC1 10 with no detectable mRNA in tumor and/or normal tissue 
Values were calculated against background; mRNA detection was categorized by GCOS software as present (P), marginally present (M), or 
absent (A). Gene symbol and title can be found at www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=gene. N/A, not applicable. 



CLIC5 


Chloride intracellular channel 5 


NLLPKLHW 


3.00 


-2.75 


A 


P 


RCC099 


ACVRL1 


Activin A receptor type ll-like 1 


SPRKGLLML 


-0.65 


3.01 


P 


A 


RCC100 


EHD2 


EH-domain-containing 2 


ALASHLIEA 


1.07 


3.32 


P 


A 


RCC099 


EHD2 
LSP1 


EH-domain-containing 2 


ALASHLIEA 


1.82 


3.41 


P 


A 


RCC100 


Lymphocyte-specific protein 1 


KLIDRTESL 




1.90 


P 


A 


RCC110 


PLVAP 


Plasmalemma vesicle-associated protein 


KVKTLEVEI 


o!50 


0.52 


P 


A 


RCC099 


SRXN1 


Sulfiredoxin 1 homolog 


TLSDLRVYL 


0.43 


0.32 


P 


M 


RCC099 


WDR78 


WD repeat domain 78 


TSWYDVAW 


0.32 


0.52 


P 


A 


RCC099 


ATOX1 


ATX1 antioxidant protein 1 homolog 


RVLNKLGGVK 


-1.01 


N/A 


A 


A 


RCC100 


CCDC21 


Coiled-coil domain-containing 21 


RLQMEQMQL 


-1.14 


N/A 


A 


A 


RCC100 
RCC099 


CYHR1 


Cysteine/histidine-rich 1 


HLGPEGRSV 


-1.21 


N/A 


M 


A 


DHX38 


DEAH box polypeptide 38 


VLFGLLREV 


0.81 


N/A 


A 


A 


RCC099 


FLJ32206 


Hypothetical protein FU32206 


GSHFISHLS 




N/A 


A 


A 


RCC099 


GBP4 


Guanylate-binding protein 4 


KRLGTLWTY 


-0.60 


N/A 


A 


A 


RCC099 


KIAA1305 


KIAA1305 


TLADIIARL 


-0.10 


N/A 


A 


A 


RCC099 


RASL11A 


RAS-like, family 11, member A 


YLLPKDIKL 


-0.37 


N/A 


A 


A 


RCC099 



histology and the broadest possible overlap of the MHC alleles 
between the different patients. Tumor samples were taken from 
its central area; control samples were from renal cortex as- 
sessed as "normal" and as distant as possible from the tumor. 

Individual Variances in Antigen Processing-ln the three 
RCC tissue pairs analyzed in this study, each sample displays 
to some extent individual features regarding antigen process- 
ing. Apart from HLA mismatches, the activity of enzymes 
involved in antigen processing also varied. We found, for 



example, varying mRNA levels of immunoproteasome sub- 
units in the different samples (Table II). All individual variances 
in antigen processing summarize at the level of the MHC 
ligandome. For example, we detected ADFP-derived peptides 
ranging in their presentation from exclusive presentation on 
tumor tissue to equal presentation on tumor and normal tis- 
sue. These individual variances in MHC:peptide levels are one 
of the reasons why for example tumor vaccination studies can 
vary extremely in their outcome. 
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Four Major Hypotheses -The HLA ligandome is more com- 
plex in its variations than mRNA, but it is of far greater immu- 
nological relevance; it reflects a particularly complex aspect of 
systems biology that is exemplified here by a naturally occur- 
ring change in tissue (tumor genesis). The long and complex 
path from mRNA transcription to HLA presentation is reflected 
in our observation of a generally weak mRNA-HLA-ligand cor- 
relation. Four hypotheses can thus be deduced from Fig. 4. 

First, HLA ligand generation implies that one mRNA species 
is able to account for several different peptides that differ in 
their presentation ratios and that are potentially presented by 
different HLA allotypes (Table IV). 

Second, for more than 75% of all identified HLA ligands, we 
could detect no major differences in the HLA peptide reper- 
toire (39) (Fig. 4, areas IV-VI), although immunoproteasomes 
were up-regulated on a transcriptional level. Assuming tran- 
scription of immunoproteasomal mRNA and activity of the 
immunoproteasomal subunits, this would limit the implication 
of the proteasome composition on the HLA ligandome in our 
ex vivo setting (36). 

Third, several HLA ligands could be identified although no 
corresponding mRNA was detectable. In these cases HLA 
ligands might be derived either from long lived proteins or 
from mutated proteins for which the correct peptide-mRNA 
allocation was not possible. 

Fourth, strong differences in HLA presentation without 
mRNA changes (Fig. 4, areas II and VIII) are due to altered 
peptide generation that can be triggered by several scenarios 

(40) . Higher translation efficiency leads to a larger number of 
DRiPs, which are known to be a major source for HLA ligands 

(41) . The amount of DRiPs from certain source mRNAs might 
also be elevated due to mutations in some of the mRNAs, 
distorting transcription and thus increasing DRiP formation 
(17). On the other hand, enhanced protein turnover increases 
the total amount of substrates for proteases such as the 
proteasome. Moreover proteolytic activity in the cytosol or ER 
can vary for some peptides when tumor and normal tissue are 
compared. 

Context and Limitations of the Chosen Experimental Set- 
ting—The exact molecular mechanism responsible for alter- 
ations in the level of our ex vivo isolated MHC-peptide com- 
plex levels could not be deduced. This is due to the fact that 
these alterations are most probably not caused by one single 
change in the peptide generation pathway or in the antigen- 
processing machinery but are rather complex, multifactorial 
changes established during tumor genesis. Defining molecu- 
lar mechanisms in a complex interaction network such as the 
antigen-processing machinery is in many cases not possible. 
In a recent study by Milner er al. (38), the turnover kinetics of 
MHC peptides was determined in a cell culture setting. Not 
even in this plain setting could the "simple" discrimination 
between short lived proteins and DRiPs clearly be made. In 
these experiments the only strong indication for MHC-peptide 
complexes derived from DRiPs was a biphasic turnover kinet- 



ics in pulse-chase experiments using incorporation of iso- 
tope-labeled amino acids. 

In 1999 and 2003 the first quantitative experiments were 
carried out to compare the yeast transcriptome with its pro- 
teome (2, 4). Here the authors showed only a weak correlation 
between transcriptome and proteome (r = 0.45). In our study, 
the correlation between transcriptome and HLA ligandome 
was even weaker (r = 0.32), which might reflect that the 
generation of HLA ligands is a pathway that begins only after 
mRNA transcription. At the present time, large scale data 
regarding the extent of correlation between proteome and 
HLA ligandome are not available. Taking into consideration 
that a major portion of the HLA ligandome is supposed to be 
derived from the transient proteome, that is short lived pro- 
teins and DRiPs, one could expect an even weaker correlation 
between the HLA ligandome and the well quantifiable perma- 
nent proteome. 

Summary and Otvf/oo/c-mRNA-based identification of tu- 
mor-associated antigens is well established and widely used. 
As we could show in this study that mRNA levels alone do not 
adequately reflect cellular reality at the HLA level, our results 
have a direct impact on T-cell-based immunotherapy. We 
conclude that for the rational selection of appropriate tumor- 
specific T-cell targets, quantitative HLA ligand analysis con- 
tributes greatly to the benefits of quantitative transcriptome 
analysis; a combination of both strategies enables the identi- 
fication of new target candidates and helps to avoid false 
positive targets, in summary, comparing changes in the tran- 
scriptome to those in the HLA ligandome of renal tumor 
versus normal tissue without further reflection is like compar- 
ing apples to oranges. Conveying this finding to RNA-based 
vaccination, which was shown to be a potent vaccine also 
against cancer (42), one cannot take for granted that every 
mRNA leads to actual peptide presentation on HLA mole- 
cules. Thus for each mRNA that is designed for vaccination, 
intensive individual in vitro tests are advisable. Our findings 
provide insight into tumor-associated changes in translation 
efficiency and protein turnover and provide immunologically 
relevant tumor information to a depth that could not be 
achieved before with hitherto existing standard techniques. 
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Currently, no method allows direct and quantitative 
comparison of MHC-presented peptides in pairs of samples, 
such as transfected and untransfected, tumorous and normal 
or infected and uninfected tissues or cell lines. Here we 
introduce two approaches that use isotopically labeled 
reagents to quantify by mass spectrometry the ratio of peptides 
from each source. The first method involves acetylation 1 and is 
both fast and simple. However, higher peptide recoveries and 
a finer sensitivity are achieved by the second method, which 
combines guanidination 2 and nicotinylation 3 , because the 
charge state of peptides can be maintained. Using differential 
acetylation, we identified a beta catenin-derived peptide in 
solid colon carcinoma overpresented on human leucocyte 
antigen-A (HLA-A)*6801. Guanidination/nicotinylation was 
applied to keratin 18-transfected cells and resulted in the 
characterization of the peptide RLASYLDRV (HLA-A*0201), 
exclusively presented on the transfectant. Thus, we 
demonstrate methods that enable a pairwise quantitative 
comparison leading to the identification of overpresented 
MHC ligands. 

The general strategy for differential MHC quantification is shown in 
Figure 1. Peptides from two different sources are derivatized either by 
an isotopically light reagent (hydrogen-containing) or an isotopically 
heavy reagent (deuterium-containing) to achieve different masses but 
the same physico-chemical behavior for the same peptide sequences. 
The two peptide derivatives are combined and separated by chro- 
matography (either offline high performance liquid chromatography 
(HPLC) or online liquid chromatography mass spectrometry 
(LCMS)). Electrospray ionization mass spectrometry (ESI-MS) is used 
to achieve greater accuracy in differential ligand quantification. 
Sequence information is then determined by tandem mass spectrome- 
try (MS/MS) and interpretation of fragment spectra using computer- 
assisted database searching tools. 

First, we established an acetylation method using a synthetic peptide 
mixture of nonamers and decamers corresponding to normal 
MHC class I ligands. The modification procedure was optimized to 
obtain complete acetylation at the N terminus. Within 15 min, we 
obtained peptides fully acetylated at the N terminus, although partial 



derivatization of the e-amino group in lysine-containing peptides also 
occurred (data not shown). 

We demonstrated the applicability of this method for differential 
acetylation of natural MHC class I ligands using an HPLC-fraction of 
HLA class I eluted peptides from MGAR cells. The peptide fraction 
was divided into two parts, one acetylated with 'Hg-acetic anhydride, 
the other with 2 D 6 -acetic anhydride, and after equimolar mixing the 
ratios of 'Hj-derivatized to 2 D 3 -derivatized peptides were determined. 
Figure 2 shows an example of the 1 H 3 -acetylated and 2 D 3 -acetylated 
variants of the peptide EVNGLISMY (MW 1040.5 Da) derived from 
the U5 small nuclear ribonucleoprotein-specific protein with a 
2 D 3 /'H 3 ratio of 1.0. For 15 additional peptides, detected as singly or 
doubly charged ions, the mean signal ratio was 1.01 and s.d. = 0.13 
(see Supplementary Table 1 online). The acetylation procedure not 
only allows differential quantification of MHC ligands from two dif- 
ferent sources but also greatly simplifies interpretation of MS/MS 
spectra of '^-acetylated and 2 D 3 -acetylated peptide species. Owing to 
the 3-Da shift of the b-series ions, N- and C-terminal fragments can be 
easily distinguished when both species are present 4 (Fig. 2c). Unlike 
tryptic fragments labeled by isotope-coded affinity tag (ICAT), most 
MHC ligands do not carry a C-terminal lysine or arginine, and there- 
fore the loss in sensitivity is quite severe for doubly charged ions 
because the acetyl group impedes the formation of a positively charged 
N terminus (signal displacement to singly charged ions occurs; data 
not shown). Usually, sequences occurring as singly charged ions yield 
signals of much lower intensity than those occurring as doubly 
charged ions (see Supplementary Table 1 online). Additionally, a sig- 
nal distribution due to partial acetylation of the s-amino group of 
lysines was observed in both 'H 3 -acetylated and 2 D 3 -acetylated pep- 
tides, resulting in lower intensity of a given signal; however, this did 
not affect quantitative comparison because both species are affected to 
the same degree under stringently controlled experimental conditions. 
Nevertheless, this method is quick, easy to perform and reliable. 

Using this differential acetylation approach, we modified peptide 
pools extracted from HLA class I molecules of solid colon carcinoma 
and regular colon tissue samples of comparable size by 2 D 3 -acetylation 
and 'H 3 -acetylation of N termini, respectively, and then combined the 
samples. After peptide separation by microbore-HPLC, nineteen natu- 
ral HLA ligands were identified by nanospray ESI-MS/MS analyses. 
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Seventeen peptides could be differentially quantified (see Supp- 
lementary Table 2 online). Most of these peptides were presented in 
comparable amounts in both tissues, yielding ratios from 1.07 to 2.42. 
The overall peptide presentation in tumor tissue was 1.7-fold higher 
than in normal tissue (mean value 1.71, Fig. 3a (left panel)). Two 
peptides showed overpresentation in tumor tissue (one example 
shown in Fig. 3a (central panel)) and one peptide was underpresented 
(Fig. 3a (right panel)). MS/MS spectra of both 1 H 3 -acetylated and 
2 D 3 -acetylated DAAHPTNVQR are shown in Figure 3c. Applying 
the student f-test to the quantified ligands, only the three over- or 
underpresented ligands are not contained within a 99.99% confidence 



Figure 1 Strategy for differential quantification of MHC-eluted peptides. 
MHC ligands derived from two different sources {e.g., tumor and normal 
tissue) are N-terminally derivatized either by a 1 H X - or 2 D„ -reagent 
(reagents: acetic anhydride; Nic-NHS ester) and combined. After reduction 
of the peptide complexity by HPLC, peptides are quantified by ESI-MS 
analysis according to their peak areas. A pair of derivatized peptides 
(H, 'Hx-derivatization; 0, 2 D„-derivatization) is physico-chemically identical 
and easily detectable because it essentially coelutes in chromatographic 
systems. Furthermore, there is a constant mass difference measured in 
the mass spectrometric scans. This difference depends on the number of 
stable isotopes in the derivative ([M+H1+: acetylation: 3 Da ; nicotinylation: 
4 Da; [M+2H] 2+ : acetylation: 1.5 Da; nicotinylation: 2 Da). Sequence 
identification of a ligand is revealed by ESI-MS/MS analysis and computer- 
assisted database search of the spectrum recorded. 



interval of 0.87 to 2.56. The two overpresented HLA ligands are 
derived from ribosomal protein L24 and beta-catenin. Although little 
data exist on the potential tumor association of ribosomal protein L24, 
it has been reported that ribosomal proteins might play a role in car- 
cinogenesis, for example L15 in esophageal carcinoma 5 . Beta-catenin 
functions as an oncogene by transactivating other oncogenes, includ- 
ing MMP-7, a metalloproteinase involved in metastasis development 6 . 
It also acts as a transcriptional coactivator in Wnt signaling, and its 
nuclear accumulation is critical for activation of the Wnt trans- 
criptional response 7-9 . Genetic manipulations that perturb nuclear 
localization of beta-catenin lead to a loss of proliferative capacity in 
colon cancer cells 10 . A number of natural HLA ligands from beta- 
catenin are listed in the SYFPEITHI database 11 (http://www. 
syfpeithi.de/). A mutated beta-catenin peptide has been defined to 
function as a target of HLA-A*24-restricted, melanoma-specific, 
tumor-infiltrating lymphocytes 12 . 

Because the acetylation method results in a substantial loss of sensi- 
tivity in peptides without C-terminal charge, we used a combination 
of two derivatization methods. For the first step, lysine side chains 
were completely guanidinated by O-methyl isourea hemisulfate 2 for a 
uniform derivatization of all s-amino groups without modification of 
the N terminus (data not shown) 13 . In a second modification step, all 
N termini were nicotinylated with nicotinoyloxy succinimide (Nic- 
NHS) reagent to maintain the charge state of the peptides 3 . This reac- 
tion was done on a Cj 8 microcolumn, which also allowed quick 
desalting before mass spectrometric analysis. Four synthetic peptides 
were used to optimize the combination of the derivatization methods 
on MHC-presented peptides. Nicotinylation of tyrosines was the only 
side reaction observed but was reversed immediately by hydroxy- 
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Figure 2 MS spectra before and after acetylation. 
(a,b) Mass spectra of nonmodified (a) and 
^-^Ds-acetylated (b) HLA peptides eluted 
from MGAR cells. A peptide fraction was divided 
into two equal parts before differential acetylation 
and afterwards combined for nanospray ESI-MS 
analysis, (c) Fragment spectra of [M+H] + of three 
species of the peptide EVNGLISMY: nonmodified, 
^-acetylated and 2 D 3 -acetylated. (d) Nomen- 
clature of fragment ions. 
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Rgune 3 Differential quantification of HLA ligands derived from two different sources, (a) Mass 
spectra of 2 D 3 -acetylated HLA peptides from colon carcinoma compared to ^-acetylated HLA 
peptides from regular colon tissue. Example of a peptide with average 203/% ratio (left); peptide 
overpresented in colon carcinoma (middle); peptide underpresented in colon carcinoma (right), 
(b) Mass spectra of guanidinated ^-nicotinylated HLA peptides from Awells cells compared to 
guanidinated 2 D 4 -nicotinylated HLA peptides from keratin 18-transfected Awells cells. Peptide 
equally presented in both cell types (left); peptide overpresented in keratin 18 transfectant (middle); 
peptide with the lowest ^/^ ratio determined in this experiment (right), (c) MS/MS spectra of !h 3 - 
acetylated (upper) and 2 D 3 -acetylated (lower) DAAHPTNVQR. (d) MS/MS spectra of 2 D 4 -nicotinylated 
synthetic RLASYLDRV (upper) in comparison to the 2 D 4 -nicotinylated natural sequence (lower). 



lamine treatment (data not shown). The strong impact of acetylation 
or guanidination/nicotinylation on peptide ionization can be seen in 
Figure 4. As indicated by the results for the three lysine-containing 
peptides, substantial loss of peptide recovery after acetylation can be 
partially compensated for by guanidination in lysine-containing pep- 
tides. The largest effect on ionization is caused by N-terminal nicotiny- 
lation, which leads to peptide recoveries similar to those of 
nonmodified peptides (Fig. 4). To determine whether guanidina- 
tion/nicotinylation is generally applicable to quantification of HLA- 
presented peptides, we used a peptide pool obtained after 
immunoprecipitation of HLA-A*0201 molecules of C1R-A2 cells. 

After guanidination of e-amines of lysine residues, this peptide pool The acetylation procedure is simple and quick compared withderiva- 
was divided into two parts, modified separately by 2 D 4 - or 'H 4 - tization by guanidination/nicotinylation. However, successful appli 
nicotinylation, respectively, and combined again. After peptide separa- tion of N-terminal acetylation is restricted to MHC-eluted peptides 
tion by microbore HPLC, 19 natural ligands were identified in the containing basic amino acids at the C-terminal anchor to maintain 
HPLC fractions by nanospray ESI-MS/MS. The mean value of signal analytical sensitivity during MS analysis. The analysis of guanidinated 
ratio for 2 D 4 Nic-^Hj Nic-peptide was 1 .00, s.d. = 0. 12, a value slightly and nicotinylated peptides is as sensitive as the analysis of nonderiva- 
better than that achieved with acetylation (data not shown). tized peptides. Thus, the two methods can be used alternatively 

To identify ligands presented exclusively by certain cells, we used according to the sample characteristics for relative quantification of 



first experiment by online-microcapillary 
LCMS. In a second MS/MS experiment, 27 
peptide sequences were identified (see 
Supplementary Table 3 online). All peptides 
were presented within the confidence interval 
on both nontransfected and transfected cells 
(Fig. 3b (left panel)), with the exception of 
one peptide with a MW of 1091.6 Da (Fig. 3b 
(central panel)). The sequence was identified 
as RLASYLDRV, which is a peptide derived 
from keratin 18. MS/MS spectra of RLASYL- 
DRV and its synthetic counterpart are shown 
in Figure 3d. There was no 'H 4 signal 
detectable in the background, indicating 
exclusive presentation of the keratin 
18-derived ligand in the keratin-transfected 
cells. For this reason an exact ratio could not 
be determined; however, the signal of the 
deuterated peptide was at least six times 
greater than the background value. The pep- 
tide with the lowest ratio is shown in Figure 
3b (right panel); however it is not signifi- 
cantly underpresented on the transfected 
Awells cells. It is contained within the confi- 
dence interval 0.64 to 2.28 obtained using the 
student f-test. 

Both the acetylation and guanidination/ 
nicotinylation methods described above are 
aimed at quantitative investigation of MHC 
ligand presentation. Since early, pioneering 
work 18 , there have been many attempts to 
detect and quantify differences in MHC-presented peptides of sample 
pairs by various mass spectrometric methods. For example, after using 
proteasome inhibitors, the pattern of HLA ligands has been compared 
to native cells 19 , and upregulation of abundant self-peptides after viral 
infection has been investigated by HPLC-MS 20 . Matrix-assisted laser 
desorption ionization/time-of-flight technology has been used to ana- 
lyze quantitative differences between peptide pools presented by 
closely related HLA allelic products 21 . None of the reported strategies, 
however, enables the fast detection of exact quantitative differences oi 
larger numbers of sequences that were previously unknown. We 
pose isotope labeling as the method of choice for such investigations. 



guanidination and 2 D 4 -/'H 4 -nicotinylation for the identification and 
quantification of HLA-ligand pools from the cell line Awells and its 
keratin 18 transfectant. Keratins have been suggested as markers that 
distinguish normal and tumor-derived cells 14 . More specifically, ker- 
atin 18 can be tumor-associated and overexpressed in carcinomas 15,16 . 
Therefore, our interest has been directed towards this protein 17 , and 
transfectants were used to identify new HLA ligands of this protein. 



MHC ligands from two different sources. 

There is a chromatographic isotope effect when using deuterium for 
stable isotope labeling: early in the elution of a peak, the deuterated 
species is enriched, whereas at the tailing edge of a peak the protonated 
peptide variant prevails. This effect increases with the number of deu- 
terium atoms contained in the reagent 22 . Thus, acetylated (D 3 ) or 
nicotinylated (D 4 ) peptides exhibit a higher degree of coelution than 



Two peptide pools extracted from HLA class I molecules of peptides containing the classical ICAT reagent with eight deuterium 

Awells-keratin 1 8 and Awells were derivatized by guanidination and atoms. A second generation of 13 C ICAT reagents has been developed 

by N-terminal 2 D 4 -nicotinylation and 1 H 4 -nicotinylation, respec- to optimize the coelution 23 . Unfortunately, reagents for the synthesis 

tively. The labeled peptide pools were combined and quantified in a of a 13 C-containing Nic-NHS reagent, which would likewise not 
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D Nonmodified ■ Aeetylatton 

■ Quanidination + acetylation ■ Guanidinstion + nicotinylation 

Figure 4 Recoveries of different peptide derivatives. The four synthetic 
peptides (AETSYVKVL, KLSLGLPGL, SLGLQLAKV and VLDPRGIYL) present 
in equimolar concentrations were either acetylated, acetylated/guanidinated 
or guanidinated/nicotinylated. After complete derivatization, each peptide 
mixture was combined with the initial nonmodified peptide mixture (thus 
expecting equimolar yields of both species) and quantified by nanospray-ESI 
MS analysis. 



exhibit a chromatographic isotope effect, are not yet commercially 
available. 

Our results show that by pair wise quantitative analysis of MHC lig- 
ands, we can detect potential tumor-associated (that is, overexpressed) 
antigens by the greater abundance of their HLA ligands and potential 
tumor-specific (that is, exclusively expressed) antigens by the existence 
of single peptide species. The identification of new tumor-associated 
antigens is usually done with quantitative screening methods such as 
gene expression analysis or proteome analysis. It is a long way, 
however, from mRNA or protein expression levels to MHC ligands, 
involving the specificity steps of proteasomal processing, transporter- 
associated-with-antigen-processing (TAP) transport and MHC bind- 
ing. Future studies of the correlation between differential quan- 
tification of MHC ligands using our approach and gene expression 
profiling or other methods would be of interest. 

METHODS 

Materials. Reagents for synthesis and derivatization, such as acetic- 'H 6 - 
anhydride, O-methyl isourea hemisulfate, nicotinic acid ethyl ester, N-hydroxy 
succinimide, EDC Ar-(3-dimethylaminopropyl)-N'-ethyl-carbodiimide hydro- 
chloride and hydroxyl amine, were purchased from Fluka. Acetic- 2 D 6 - 
anhydride was purchased from Sigma-Aldrich, 2 D 4 -nicotinic acid ethyl ester 
was purchased from LGC Promochem. The HPLC reagents, trifluoroacetic 
acid, acetonitrile, formic acid and HPLC water, were purchased from Merck. 

Synthesis of the Hf'H^DJ nicotinoyloxy)succinimide (>H 4 / J D 4 -Nic-NHS) 
esters. 'D 4 -nicotinic acid ethyl ester (5 g, 32.2 mmol) was dissolved in 
ethanol/water (1:1, 32 ml). Sodium hydroxide (1.53 g, 38.1 mmol) was added 
and the solution refluxed over night. After removal of ethanol by evaporation 
the solution was diluted with water (40 ml) and extracted with diethylether 
(32 ml). The aqueous phase was neutralized with concentrated HC1 and pH was 
adjusted to 4-4.5. The resulting precipitate was filtered and washed with cold 
water and ether. Aqueous washings and mother lye were combined, adjusted to 
pH 2.8 and evaporated to dryness. The residue was redissolved in 1 N HCl 
(as little as possible) and adjusted to pH 3.5 with saturated NaOH. The result- 
ing precipitate was filtered and washed with cold water and ether. Combination 
of the precipitates yielded ^-nicotinic acid as a white powder (3. 1 1 g, 76%). 
Thus obtained, ^-nicotinic acid (3.11 g, 24.5 mmol) was suspended in 
dry CH 2 C1 2 (150 ml). JV-(3-dimethyIaminopropyl)-N'-ethyI-carbodiimid 
hydrochloride (EDC) (5.18 g, 27 mmol) and N-hydroxy succinimide (3.53 g, 
30.63 mmol) were added and the resulting solution stirred overnight at 
22 C. The reaction mixture was extracted with saturated sodium bicarbonate 



solution (75 ml) and water (75 ml) and the organic phase was dried (NaS0 4 ). 
Evaporation of the solvent in vacuo yielded the 1 H 4 / 2 D 4 -Nic-NHS ester (5.06 g, 
92%). Purity and structure were confirmed by 'H- and 13 C-NMR. 

Peptides. Synthetic peptides were synthesized in an automated peptide synthe- 
sizer 432A (Applied Biosystems) following the 9-fluorenylmethyl-oxocar- 
bony/fert-butyl (Fmoc/tBu) strategy as described 24 . 

Elution of HLA-presented peptides. HLA peptides were obtained by immune 
precipitation of HLA molecules from cell lines (HLA-defined human B-lym- 
phoblastoid lines C1R-A2, International Histocompatibility Workshop (IHW) 
No. 9208; MGAR, IHW No. 9014; Awells, IHW No. 9090) or solid tissues 
according to a slightly modified protocol 25 - 26 using the HLA-A-,-B-,-C-specific 
antibody W6/32, protein A-sepharose, acid treatment and ultrafiltration. 

Transfection. The cDNA of keratin 18 was subcloned using TOPO TA cloning 
(Invitrogen) and inserted into the EcoRI and Nofl sites of pcDNA3-Ii, in frame 
with the Ii sequence. Stable transfectants were generated by electroporation of 
Awells cells, followed by cloning using the limiting dilution method. 

Acetylation of peptides. We added 10 pi 'H 6 -acetic anhydride or 2 D 6 -acetic 
anhydride (50% (vol/vol) in methanol) to 100 pi peptide mixture (200 pmol to 
2 nmol) in 50% methanol/water (vol/vol). The mixture was allowed to react at 
22 Cforl5min.The reaction was stopped by addition of 1 . 1 pi of formic acid 
and equal aliquots of both samples were mixed. 

Guanidination of peptides. Peptide mixtures (200 pmol to 2 nmol) in citrate 
buffer (50 mM; pH 3) containing 0.25% trifluoroacetic acid (TFA) (vol/vol) 
were adjusted to a pH of 10.5 with 200 pi sodium hydroxide ( 10 M). After addi- 
tion of 1 ml freshly prepared O-methyl isourea hemisulfate solution (2.5 M in 
water) the mixture was incubated for 10 min at 65 C in a water bath. The reac- 
tion was stopped by adding 200 pi formic acid. 

Nicotinylation of guanidinated peptides. Guanidinated peptides were loaded 
on a reversed phase C-18 microcolumn (Agilent Technologies hydrophobic 
202-XGSXB) and washed with 0.5 ml water. The bound peptides were then 
modified by slowly passing 1 ml of freshly prepared 'H 4 - or 2 D 4 -Nic-NHS ester 
(sodium phosphate buffer 50 mM; pH 8.5) through the column. The procedure 
was done twice before 1 ml of hydroxylamine was run through the column to 
remove tyrosine modifications. Finally, after washing with water, the modified 
peptides were eluted using 100 pi of 50% acetonitrile/water (vol/vol). 

Offline HPLC separation. Modified peptide pools were mixed in equimolar 
amounts and the volume reduced to approximately 100 pi by vacuum centrifu- 
gation. The mixture was diluted with 400 pi water containing 0.08% TFA 
(vol/vol) before injecting the sample into a microbore SMART-HPLC system 
(Amersham). Labeled peptide mixtures were separated by gradient elution 
from a reversed-phase column (pRP SC C2/C18-column, 100 2.1 mm, 
Pharmacia). The samples were run through a binary gradient of 10-60% B 
within 100 min, applying a flow rate of 150 pl/min. Solvent A was 0.1% 
TFA/water (vol/vol). Solvent B contained 0.08% TFA in 80% acetonitrile/water 
(vol/vol). The peptides were fractionated in 1 50 pi aliquots. Before MS analysis 
each fraction was dried completely and resuspended in 50% 
methanol/water/0.1% formic acid (vol/vol). 

Microcapfflary ICMS. Synthetic and natural peptide mixtures were analyzed 
by a reversed phase Ultimate HPLC system (Dionex), coupled to a hybrid 
quadropole orthogonal acceleration time-of-flight MS/MS (Q-TOF, 
Micromass) 27 equipped with a micro-ESI source. Samples were loaded onto a 
C 18 precolumn for concentration and desalting. After loading, the precolumn 
was placed in line for separation by a fused-silica microcapillary column 
(75 pm internal diameter 250 mm) packed with 5 pm C, 8 reversed-phase 
material (Dionex). Solvent A was 4 mM ammonium acetate/water. Solvent B 
was 2 mM ammonium acetate in 80% acetonitrile/water. Both solvents were 
adjusted to pH 3.0 with formic acid. Samples were run through a binary gradi- 
ent of 15-60% B within 120 min, applying a flow rate of 200 pl/min reduced to 
~300 nl/min by the Ultimate split-system. A gold-coated glass capillary 
(PicoTip, New Objective) was used for introduction into the micro-ESI source. 
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The integration time for the TOF analyzer was 1 s with an interscan delay of 
0.1 s. The ratios of deuterated and nondeuterated peptide species were calcu- 
lated from the relative peak height. 

For online microcapillary HPLC MS/MS experiments, the integration time 
for the TOF analyzer was 4 s with an interscan delay of 0.1 s. Fragmentation of 
parent ions [M+H] + or [M+H] 2+ was achieved by automatic switching 
between MS and MS/MS. Fragment spectra were analyzed manually and data- 
base searches (National Center for Biotechnology Information, Expressed 
Sequence Tag) were carried out using Multiple Alignment System for Protein 
Sequences Based on Three-way Dynamic Programming (MASCOT). 

Nanoflow ESI-MS analysis and quantification. For sample introduction metal- 
coated glass capillaries (Proxeon) were used, working at flow rates of about 
20 to 50 nl/min. The ratios of deuterated and nondeuterated peptide species 
were calculated from the relative peak heights as well as from peak area. For 
Nanoflow ESI MS/MS experiments, peptide ions were selected and fragmenta- 
tion of [M+H] + and [M+2H] 2+ was carried out using collision energies of 
30-60 eV and 20-30 eV, respectively. The integration time for the TOF analyzer 
was 1 s with an interscan delay of 0.1 s. 

Note: Supplementary information is available on the Nature Biotechnology website. 
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